
Braz J Med Biol Res, November 2012, Volume 45(11) 1080-1085 



doi: 1 0. 1 590/SO 1 00-879X20 12007500 1 3 1 



Association of sulfur dioxide exposure with circulatory system 
deaths in a medium-sized city in Brazil 

C.T. Amancio and L.F.C. Nascimento 



The Brazilian Journal of Medical and Biological Research is partially financed by 



Ministerio 
da Ciencia eTecnologia 




Ministerio 
da Educagao 



u i _ 

IM PAIS DE TODOS 

GOVERNO FEDERAL 



CCSP 



Campus 

Ribeirao Preto Faculdade de Medicina 
de Ribeirao Preto 



Institutional Sponsors 



SHIMADZU 





UNICAMP 



Associagao Explore High - Performance MS 
Titntivo Orbitrap Technology 

aPesquisa I" Proteomics & Metabolomics 



analitica Thermo 

analiticaweb.com.br scientific 



EY-Ni 



All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License 



Brazilian Journal of Medical and Biological Research (2012) 45: 1080-1085 
ISSN 1414-431X 



Association of sulfur dioxide exposure 
with circulatory system deaths in a 
nnediunn-sized city in Brazil 

C.T. Amancio'' and L.F.C. Nascimento'''^ 

'I Departamento de Medicina, Universidade de Taubate, Taubate, SP, Brasil 
^Programa de Pos-Graduagao em Ciencias Ambientais, Universidade de Taubate, Taubate, SP, Brasil 



Abstract 

There is a demonstrable association between exposure to air pollutants and deaths due to cardiovascular diseases. The objec- 
tive of this study was to estimate the effects of exposure to sulfur dioxide on mortality due to circulatory diseases in individuals 
50 years of age or older residing in Sao Jose dos Campos, SP. This was a time-series ecological study for the years 2003 to 
2007 using information on deaths due to circulatory disease obtained from Datasus reports. Data on daily levels of pollutants, 
particulate matter, sulfur dioxide (SO2), ozone, temperature, and humidity were obtained from the Sao Paulo State Environ- 
mental Agency. Moving average models for 2 to 7 days were calculated by Poisson regression using the R software. Exposure 
to SO2 was analyzed using a unipollutant, bipollutant or multipollutant model adjusted for mean temperature and humidity. The 
relative risks with 95%CI were obtained and the percent decrease in risk was calculated. There were 1928 deaths with a daily 
mean (± SD) of 1 .05 ± 1 .03 (range: 0-6). Exposure to SO2 was significantly associated with mortality due to circulatory disease: 
RR = 1 .04 (95%CI = 1 .01 to 1 .06) in the 7-day moving average, after adjusting for ozone. There was an 8.5% decrease in risk 
in the multipollutant model, proportional to a decrease of SO2 concentrations. The results of this study suggest that residents 
of medium-sized Brazilian cities with characteristics similar to those of Sao Jose dos Campos probably have health problems 
due to exposure to air pollutants. 

Key words: Air pollution; Sulfur dioxide; Cardiovascular diseases; Mortality; Air pollutants; Stroke 



Introduction 

Circulatory diseases are the leading cause of death in 
Brazil and are responsible for rennarkable financial and social 
costs, as observed in the State of Sao Paulo, where R$500 
nnillion (US$1 = ~R$1.70) were spent in hospitalization of 
individuals aged 50 years or nnore in 2010 due to these 
causes (1 ). Nunnerous risk factors for circulatory diseases 
have been well elucidated, such as age, diabetes, dyslipi- 
dennia, snnoking, alcoholisnn and, nnore recently, exposure 
to environnnental factors such as air pollution (2,3). 

Environnnental pollution is an innportant topic today be- 
cause it is a preventable public health problenn. Its principal 
sources of ennissions are the industries and connbustion 
nnotor vehicles and the nnain air pollutants are particulate 
nnatter (PM10), sulfur dioxide (SO2), ozone (O3), carbon 
nnonoxide (CO), and nitrogen oxides (NOx) (4). It has been 
shown that, even when within acceptable standards accord- 
ing to environnnental legislation, exposure to pollutants is a 



risk factor for both cardiovascular and respiratory nnorbidity 
and nnortality (5-8). 

It is known that the innpact of pollutants on vascular 
diseases is greater annong adult and the elderly, while 
respiratory diseases are nnore connnnon annong children 
and the elderly (9). 

Moreover, the adverse effects of exposure to air pollut- 
ants occur both in populations living in large cities and in 
nnediunn-sized cities (10-12). 

Study areas in Sao Paulo have shown an association 
between exposure and hospital adnnissions for acute nnyo- 
cardial infarction, showing a stronger association with SO2 
(10). A European study (APHEA-II) suggested that SO2 
can have an independent role in triggering cardiovascular 
events and hospitalizations due to these causes, coincid- 
ing significantly with increases in daily levels of SO2 on the 
sanne day and on the day before exposure (13). 
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Are the effects of exposure to SO2 the same in medium- 
sized cities? Faced with this challenge, the aim of this study 
was to estimate these effects on mortality due to circulatory 
disease in Sao Jose dos Campos, SP, Brazil. 

Material and Methods 

This was an ecological time-series study in which in- 
formation on deaths from circulatory disease was obtained 
for individuals 50 years of age or older. The study was 
conducted in Sao Jose dos Campos, a medium-sized city 
located 80 km from Sao Paulo, between two mountain 
ranges, with a population of about 700,000 inhabitants. 
The city is located at latitude 23° 11' S and longitude 45° 
53' W and its altitude is 600 m above sea level. It has about 
1 1 00 industrial establishments with emphasis on automobile 
manufacturing, aerospace and pharmaceutical industries, 
and oil refinery. The city is crossed by the Dutra Highway, 
which is the most important highway in Brazil with heavy 
car, truck and bus traffic involving about 130,000 vehicles 
per day (14). 

Data about the selected circulatory diseases culminat- 
ing in death were obtained from the Brazilian Mortality 
Information System (SIM). According to the International 
Classification of Diseases, 10th revision, they were coded 
as ICD-10 = no. 111, 115, and 120 to 125 (heart diseases), 
and 160 to 167 (stroke). The study considered the data for 
the period from January 1 , 2003 to December 31 , 2007. 

We obtained data about the daily levels of air pollut- 
ants such as SO2, PM10 and O3 from the Sao Paulo State 
Environmental Agency (CETESB), which has a monitor- 
ing station in the central region of Sao Jose dos Campos 
located 2 km north of the Dutra Highway. Data collection 
about the levels of pollutants was started at the first hour 
of the day and continued for 24 h. All data were quantified 
in ijg/m^, considering the daily average for SO2 and PM10 
and a maximum of 1 h of the day for O3. 

In Sao Jose dos Campos, the PM10 pollutant is com- 
posed of nitrate, sulfate and chloride anions and constituent 
particles of soil and soot; sulfate was predominantly found in 
the fine fraction (PM<2.5), the lowest values were detected 
for percent chloride ion, which was the predominant coarse 
fraction (PM2.5-10); for nitrate, the largest percentages were 
found in the coarse fraction (PM2.5-10), approximately two 
times higher than those found in the PM<2.5 (1 5). 

Temperature and relative humidity data were obtained 
from the Foundation for Science, Technology and Space 
Applications (FUNCATE). 

The dependent variable was death due to circulatory 
diseases; the independent variable was sulfur dioxide, 
adjusted for PM10, O3, mean temperature and humidity. 

Since the manifestations of exposure to pollutants are 
believed to present a lag effect, i.e., that an individual exposed 
to pollution today may present health problem on the same 
day or some days later, models with moving average ranging 



from 2 to 7 days after exposure to the air pollutant were con- 
structed. The moving average, in turn, is a tool used to obtain 
the cumulativeeffect of concentrations of pollutants and is more 
consistent with reality. Studies do not show consensus about 
the days included in the moving average (16). 

Statistical analysis 

Time-series analysis was used, with the unit of obser- 
vation being the day and not the individual. The analysis 
assessed the fraction of the daily variations in death counts 
that was explained by the daily variations in air pollution, 
SO2 in particular, of the preceding days, after controlling 
for the other variables that varied in time. We used gener- 
alized linear models of Poisson regression to estimate the 
association between exposure to SO2 and deaths due to 
circulatory system conditions in persons aged 50 years or 
more, using the R software (R Systems International, India). 
We chose this statistic approach because the unit of study, 
death due to circulatory disease, is a counting event. 

The analysis was based on a unipollutant model for 
SO2 adjusted for average temperature and humidity. Bi- 
pollutant models were constructed, a model with SO2 and 
PM10 and the other with SO2 and O3. The multipollutant 
model included SO2, particulate matter and ozone adjusted 
for temperature and humidity. Relative risks (RR) were 
estimated for diseases of the circulatory system, for heart 
disease and stroke. 

Thus, we obtained the coefficients (coeff) and its stan- 
dard deviations (SD) that permitted the calculation of RR and 
95% confidence intervals (95%CI), through the formulas: 
RR = exp (coeff) and 95%CI = exp [coeff ± 1 .96 (SD)]. 

We report descriptive analysis of all study variables 
and correlations among them through Pearson's correla- 
tion coefficients, for which we used the computer program 
SPSS version 15.0. The decrease in risk of death for the 
interquartile difference for SO2 was estimated using the 
formula: PD = [exp (-coeff * VIQPOL) - 1] * 100, where PD 
is the percent decrease in the risk of circulatory deaths and 
VIQPOL is the difference between the values of the first and 
third quartiles of the concentration of this pollutant. PD was 
estimated only for SO2, considering the difference between 
the values of the first and third quartiles of the concentra- 
tion of this pollutant obtained by the multipollutant model 
for circulatory diseases. The level of significance was set 
at 5% for all tests. 

Results 

A total of 10,111 deaths of individuals aged 50 years or 
more occurred during the study period from January 1 , 2003 
to December 31 , 2007; 1 928 of these deaths (19.1 %) were 
due to circulatory diseases. There were 1 021 (53.0%) deaths 
due to cardiovascular diseases and 907 (47.0%) deaths due 
to stroke. The descriptive analysis of the variables under 
study is given in Table 1 . The average particulate matter was 
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Table 1. Daily mean values of the variables studied in Sao Jose 
dos Campos, Brazil, 2003-2007. 



Table 2. Pearson's correlation matrix between all atmospheric 
variables. Sao Jose dos Campos, Brazil, 2003-2007. 





Mean ± SD 


Minimum 


Maximum 


Deaths 


1.06 ± 1.03 


0 


6 


PMio (Mg/m^) 


26.97 ± 14.8 


4 


100 


SO2 (Mg/m3) 


4.38 ± 3.52 


1 


34 


O3 (MQ/m^) 


80.00 ± 33.3 


5 


232 


Temperature (°C) 


21.73 ±3.84 


12.7 


36.2 


Humidity (%) 


76.24 ± 15.0 


33 


100.0 



PM10 = particulate matter. 

26.9 ± 14.8 ijg/m^, not exceeding the established pattern 
of 50 ijg/nn^ annual average or the average annual daily 
average of 150 |jg/nn^. Annual average daily level of SO2 
was 4.4 ± 3.5 |jg/nn^ and a nnaxinnunn of 1 h of the day for 
O3 was 232 ± 33.3 ijg/nn^, with acceptable values for the 
two pollutants having been established up to 80 and 160 
pg/nn^, respectively. This shows that SO2 did not exceed 
the acceptable value, whereas O3 exceeded this linnit on 
34 occasions for a total of 1826 days of study. 

The year-by-year variation during the study period is 
shown in Figure 1 . There was a high seasonal behavior for 
PM and SO2, but not for O3. 

The nnortality rate fronn all circulatory diseases varied 
little in ternns of total annual nunnber, with the lowest nunn- 
ber occurring in 2005 (366 deaths) and the largest in 2007 
(401 deaths). However, when the rate was analyzed on 
a nnonthly basis there was seasonal variation, with larger 
nunnbers in the colder nnonths of the year and lower nunnbers 
in warnner nnonths. 

Table 2 presents the correlation nnatrix between the 
variables of the study. We observed positive correlations 
between pollutants, whereas deaths due to circulatory 
diseases showed a positive correlation with PM10 and SO2, 
and a negative correlation with O3. The weather variables 
showed negative correlations with deaths and air pollution, 
except for tennperature, which was positively correlated 
with PM10 and O3. 

A generalized linear nnodel using Poisson regression 
analysis was used to deternnine the association between 
exposure to SO2 and deaths due to circulatory diseases and 
is expressed as RR (Figure 2A-D). The nnodel presents the 
RR and their 95%CI for SO2 in the unipollutant nnodel, in the 
bipollutant nnodel and in the nnultipollutant nnodel with PM 
and O3, adjusted for nnean tennperature and hunnidity. The 
RR for death due to circulatory diseases, heart diseases 
and stroke are shown in Table 3; it is possible to note the 
contribution of death due to stroke to nnortality due to cir- 
culatory diseases. The nnaxinnunn RR obtained was 1.04 
(95%CI = 1 .01 to 1 .06) for the bipollutant nnodel adjusted for 
O3. There was a decreased risk of death of approxinnately 
8.5% for the fifth, sixth, and seventh nnoving average, in 
the nnultipollutant nnodel for a decrease in the interquartile 





Deaths PM10 


SO2 


O3 


Temperature 


Humidity 


Deaths 


1 0.05* 


0.06* 


-0.02 


-0.04 


-0.01 


PM10 


1 


0.51* 


0.39* 


0.16* 


-0.37* 


SO2 




1 


0.26* 


-0.03 


-0.30* 


O3 






1 


0.04* 


-0.35* 


Temperature 








1 


-0.68* 


Humidity 










1 



PM10 = particulate matter. *P < 0.05. 
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Figure 1. Daily pollutant concentration during the study period. Sao 
Jose dos Campos, Brazil, 2003-2007. A, PM10; B, SO2, C, O3. 

difference (3 pg/nn^) of SO2 (Figure 2E). 

Discussion 

This is the first tinne-series study reported to estinnate 
the association between exposure to air pollutants and 
circulatory deaths in a nnediunn-size city in Brazil. 

This was an ecological study, since the unit of study 
was the population and not the individual. The choice of 
statistical analysis using a generalized linear nnodel rather 
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Figure 2. Relative risks and 95% confidence interval for deaths 
due to cardiovascular circulatory diseases by exposure to sul- 
fur dioxide: A, unipollutant model; B, bipollutant model adjusted 
for PMio; C, bipollutant model adjusted for O3; D, multipollutant 
model adjusted for PM10 and O3; E, percent decrease in risk ac- 
cording to the interquartile difference (3 |jg/nn^) of SO2 in the mul- 
tipollutant model according to moving average of two days (MA 2) 
up to 7 days (MA 7). Sao Jose dos Campos, SP, 2003-2000. 



Table 3. Relative risk (RR) and respective 95% confidence in- 
terval (95%CI) for deaths due to circulatory diseases, cardiac 
diseases and stroke in relation to exposure to sulfur dioxide in a 
multipollutant model. Sao Jose dos Campos, Brazil, 2003-2007. 



Circulatory diseases Cardiac diseases 
[RR (9S%CI)] [RR (9S%CI)] 



Stroke 
[RR (9S%CI)] 



MA 2 
MAS 
MA 4 
MAS 
MA 6 
MA 7 



1.016 (0.998-1 .OSS) 
1.019 (0.998-1.040) 
1.025 (1.003-1.048) 
1.030 (1.007-1.054) 
1.030 (1.006-1.055) 
1.030 (1.005-1.056) 



1 .011 (0.986-1 .037) 1 .022 (0.996-1 .048) 

1 .011 (0.982-1 .040) 1 .02s (0.99S-1 .OSS) 

1 .020 (0.990-1 .0S2) 1 .027 (0.99S-1 .0S9) 
1 .02s (0.993-1 .0S8) 1 .031 (0.997-1 .066) 

1.021 (0.993-1. 0S6) 1.035(1.000-1.072) 
1.019 (0.98S-1.0SS) 1.036(1.000-1.074) 



MA = moving average. P < 0.05 for the values in bold (Poisson re- 
gression). 



than a generalized additive nnodel was based on a study 
conducted by Conceigao et al. (17), which dennonstrated that 
both nnodels yield consistent results, without connpronnising 
the end result of the study. The generalized linear nnodel 
was also used successfully by Arbex et al. (18) to assess 
the association between exposure to total suspended par- 
ticles resulting fronn bionnass burning and hospitalizations 
for hypertension. 

The levels of pollutants reported in the present study dif- 
ferfronn those reported for others, such as a study conducted 
in Sao Paulo, the Brazilian nnetropolis, which dennonstrated 
24-h PMi 0 values exceeding by several tinnes the linn it value 
of 1 50 ijg/nn^ (1 9). This finding is consistent, since Sao Jose 
dos Cannpos is a nnediunn-sized city with a snnaller vehicle 
fleet (230,000 vehicles versus 5 nnillion in Sao Paulo) and 
industrial park, which are nnainly responsible for the ennission 
of PM10 (20). Also in Sao Paulo, the average SO2 value 
was 17.71 |jg/nn^, which is well above the average found 
in the present study (4.4 ijg/nn^), and O3 levels were close 
to those detected in the present study (80 ijg/m^), with an 
average of 71 .79 ijg/nn^ average (6). 

Both environnnental pollutants and the sunn of nnonthly 
deaths throughout the study period showed nnarkedly 
seasonal behavior, with increased levels of pollutants and 
nunnber of deaths in the cooler tinnes of the year (nnonths of 
May, June and July) and decreased values in the warnner 
nnonths of the year. This can be explained by the thernnal 
inversion that occurs during the cold nnonths, reducing 
the dispersion of air pollutants. Ozone, on the other hand, 
showed no seasonal pattern, since it is known that it is the 
result of chennical reactions involving sunlight, nitrogen 
oxides, and hydrocarbons derived prinnarily fronn vehicle 
ennissions, unlike the nnechanisnn of the other three pol- 
lutants. 

The nnechanisnns involved in the genesis of circulatory 
diseases by exposure to pollutants have not been well 
explored, but a possible nnechanisnn could be the increase 
in plasnna fibrinogen and inflannnnatory factors, which lead 
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to increased blood viscosity, resulting in clinical cardiovas- 
cular events (21). 

In particular, SO2 resulting from the burning of fossil 
fuels is a major component of air pollution in many parts 
of the world. In urban areas, SO2 originates from domestic 
heating, power generation through thermal power plants 
and motor vehicles. Its inhalation is responsible for adverse 
health effects due to its absorption by the mucous mem- 
branes of the respiratory tract (22). A study conducted by 
Tunnicliffe et al. (23) using electrocardiograms suggested 
that exposure to SO2 can influence the autonomic nervous 
system, contributing to the explanation of the mechanisms 
involved in the cardiovascular and bronchial constriction 
produced by this pollutant. Another study also suggested a 
direct effect on the autonomic nervous system rather than 
effects secondary to a systemic inflammatory response, 
having found that short-term exposure to SO2 causes a 
decrease in cardiac vagal control measures (24). 

In the present study, we demonstrated the association 
of exposure to SO2 and circulatory events that result in 
death. Exposure to SO2 was significantly associated with 
mortality due to circulatory diseases (maximum RR in the 
moving average of 7 days = 1 .036; 95%CI = 1 .014-1 .060), 
after adjusting for O3. Although the RR estimated in this 
study were of low magnitude, it is important to note that 
exposure to air pollutants is a frequent event, resulting in 
a marked impact on the health system and the exposed 
individual. 

In another study carried out in Sao Paulo, an increase 
of 10 ijg/m^ in SO2 concentrations led to an 11% increase 
in the risk of hospital admissions due to respiratory disease, 
a 3% increase in hospital admissions due to cardiovascular 
diseases, and a 7% increase in hospital admissions due to 
ischemic heart disease in elderly people (18). 

Along this line, Sunyer et al. (13) found in seven Euro- 
pean regions a significant increase in the daily numbers of 
all cardiovascular admissions except for stroke, particularly 
ischemic heart disease (IHD), with an increase in the levels 
of SO2 in the same day and the day before. After adjustment 
for PM10, the combination of SO2 with IHD admissions re- 
mained significant among individuals younger than 65 years, 
but not among individuals older than 65 years. An increase 
of 10 ijg/m^ of SO2 in the daily average corresponded to 
0.7% of all cardiovascular hospitalizations on the same 
day and the next. On the other hand, a study conducted in 
Finland identified an association between exposure to PM 
and death from stroke (25) 

In the present study, we chose to estimate the decrease 
in RR for a decrease in the interquartile difference in SO2 
(3 ijg/m^); under these conditions there was a decreased 



risk of death reaching a maximum value of approximately 
8.5% in the multipollutant model. This finding reinforces the 
need to implement measures to decrease the concentra- 
tions of air pollutants. 

This study may have limitations where the individual 
exposures are not considered, taking as the basis a homo- 
geneous atmosphere throughout the analysis. Therefore, 
we cannot say that an individual who died was necessarily 
exposed to higher pollution levels. The use of the general- 
ized linear model instead of a generalized additive model 
may have been another limitation, because the general 
additive model is a more parsimonious model that requires 
a smaller number of explanatory variables. Another limita- 
tion may have been a mistake in the coding of cause of 
death, besides the lack of information on comorbidities. The 
database used in this study does not identify whether the 
subject is a smoker or non-smoker. Another limitation was 
the lack of separation between deaths due to heart disease 
and deaths due to stroke. In our study, stroke made the most 
expressive contribution to the deaths recorded. 

We should be careful regarding the exclusive involve- 
ment of SO2 in the occurrence of deaths because other 
pollutants such as CO, NO and NO2 were not included 
due to the fact that they are not estimated by CETESB. It 
is important to note that this statistical approach and this 
type of study may suggest a relationship, but not neces- 
sarily a causal effect. However, it serves as a basis for the 
management of public health, social and financial aspects 
of the country. 

The results of this study suggest that residents of 
medium-sized Brazilian cities with characteristics similar 
to those of Sao Jose dos Campos, which is crossed by a 
highway with considerable truck and bus traffic, and hav- 
ing a similar industrial complex, may also be affected by 
exposure to environmental pollutants, even when they are 
present in low concentrations and within the limits accept- 
able by environmental legislation (26). 

Moreover, the death outcome examined in this study, 
although being the most serious, is just one of the outcomes 
of exposure to air pollutants, also represented by emer- 
gency consultations, hospitalizations and poorer quality 
of life. Thus, we observed an important role of SO2 when 
combined with the risk of circulatory deaths in the elderly 
in a Brazilian medium-sized city. 
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